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Abstract 

In an R-parity nonconserving super symmetric theory, the lighter stop can 
dominantly decay into b\i and br if R-parity breaking has to explain the 
neutrino mass and mixing pattern suggested by the data on atmospheric muon 
neutrinos. This should give rise to dilepton + dijet and single — lepton+ jets, 
signals identical with those of the top quark at the Fermilab Tevatron. One 
can thus constrain the stop parameter space using the current top search data, 
and similarly look for the first signals of supersymmetry at the upgraded runs 
of the Tevatron. 
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It has been repeatedly suggested that in a supersymmetric (SUSY) scenario [jlj, one of the 
two spin-zero superpartners of the top quark, popularly known as stop, could be considerably 
lighter than any other strongly interacting superparticle || . This happens due to the mixing 
between the left-and right chiral stops, which can be quite large in principle, reducing the 
mass of the lighter eigenstate. Also, negative contribution from Yukawa coupling in the 
evolution equations for scalar mass parameters can cause the stop to be lighter than the 
other squarks in a framework where all squark flavours have the same mass at a high energy 
scale. It has been shown that in the minimal supersymmetric standard model (MSSM), a 
light stop (ti) will decay dominantly through the channels t\ — > cx\ and t\ — > bxf , where Xi 
and xt are the lightest neutralino and the lighter chargino respectively. Out of these, the 
former reigns supreme || if the light stop is even lighter than xt @- Based on the ensuing 
signals, experiments at the Fermilab Tevatron have set a lower limit of about 120 — 135 GeV 
on the lighter stop || depending on the mass of the lightest neutralino. 

The observed signals can be quite different when R-parity, defined as (— \Y B+L+2S , is not 
a conserved quantity, something that is both theoretically and phenomenologically consistent 
in supersymmetric theories so long as only one of baryon number (B) and lepton number 
(L) is violated. In such cases, the lightest supersymmetric particle (LSP) can be unstable, 
thereby altering the signals of SUSY ||[7]]. If, in addition, the LSP (which is the lightest 
neutralino in most theories) is massive enough to evade searches at the large electron-positron 
collider (LEP), then the observation of SUSY at the Fermilab Tevatron may depend squarely 
on the production and decays of the light stop for whom R-parity violation opens up new 
and often dominant decay modes. The limits based on t\ — > cxi require re-examination in 
such a case. 

In this letter, we want to point out that the easiest identifiable signatures of a light stop 
(and therefore perhaps of SUSY itself) in an R-parity nonconserving framework could be 
in final states that have been, and still are, widely studied to look for the top quark at the 
Tevatron, namely, the dimuon+dijet as well as single muon+jets signals. We further argue 
that stop decays leading to such final states are imperative in R-parity violating scenarios 
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that lead to neutrino masses and mixing as required by the recent SuperKamiokande (SK) 
data on atmospheric muon neutrinos. Thus a careful reanalysis of the already existing 
Tevatron Run I results on top search may allow us to explore a large region of the light stop 
parameter space of such theories, where the search channel employed for MSSM is not going 
to be effective. In Run II one can have an even wider prospect, not only for constraining 
the SUSY parameter space but also for the possibility of top signals being actually faked 
by the stop. In addition, some new signals for such models, very similar to those of the top 
quark, are suggested here. 

Expressed in terms of the quark, lepton and Higgs superfields, the MSSM superpotential 



is 



W MSSM = ixE x E 2 + ti^UHyE] + h%Q t E x b) + h^QiHzU? (1) 

where \i is the Higgsino mass parameter and the last three terms give all the Yukawa inter- 
actions. 

The possible additions to this superpotential due to R-parity violation (through lepton 
number violation only) are given by || 

W# = tiLiH 2 + XijkLiLjEl + X ijk LiQjDl (2) 

We consider only the effects of the trilinear additional terms in the superpotential; we shall 
point out at the end of the paper that the effects we are suggesting can arise also from the 
bilinear terms e^I/ji^- 

If the SK data || on atmospheric muon neutrinos (and also similar data from the Soudan- 
II and MACRO |IJ experiments) have to be explained in terms of — v T oscillations, 



then the mass-squared splitting between the second and third lightest physical neutrino 
states will have to be Am^ ~ 5 x 1CT 3 , in addition to near-maximal mixing between 
the corresponding flavour eigenstates. R-parity violation in the form of the A-and A -type 
interactions in equation 2 can give rise to neutrino mass terms at one-loop level, the generic 
expression for them (in the flavour basis) being 
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(mt oop )y ~ —m d k m d p M SU sY^Kk P ^pk 

+ ^m l k m l p M S usY^kk P X jP k (3) 

where m d ^ denote the down-type quark (charged lepton) masses, m~, m~ are the average 
slepton and squark mass squared. Msusy^ fJ-) is the effective scale of supersymmetry 
breaking. The mass and mixing patterns suggested by the observed data can be accom- 
modated in the above scheme of mass generation if, with Msusy — mq — 300 — 500 GeV, 
A 233 — A333 ~ a few times 10~ 4 flUl - An immediate consequence of such couplings is the 



decay of the lighter stop in the channels t\ — > br + and t\ — > bfi + with comparable widths. 

There is a considerably large region of the parameter space, not yet constrained by any 
experimental data, where the stop is the second lightest supersymmetric particle next to 
the lightest neutralino. In this region, the three lowest-order decay channels available to the 
stop are ti — > CX\, h — > br + and t\ — > 6/1 + The first one, a well-studied process, is a 



consequence of neutral flavour violation in SUSY and is suppressed by the small mismatch 
between quark and squark mass matrices. The latter ones are driven by A 233 and A 333 , and 
we find that they dominate for a light stop with mass < 150 GeV so long as A' 233 , A 333 lie in 
the range specified above, in conformity with the SK data. 

We focus our attention on single-muon and dimuon final states, together with jets and 
missing energy, arising from (light) stops pair-produced in pp collisions at the Fermilab 
Tevatron, with either both decaying into br or one of them into br while the other goes 
to b/i. The same final states have been extensively analysed for the determination of the 
mass and production cross-section of the top quark. It should be noted that the r pro- 
duced together with a b quark in the two-body decay of a stop (with mass 100 — 150 GeV) 
can have sufficient pt for the resulting jets/muons/neutrinos to often pass the px or /E? 
cuts associated with the top quark signals. Thus, depending upon whether the tau decays 
hadronically or semileptonically, the stop decays can contribute to top-like signals of the 
form {%) single muon + 3jets+ /Et and (ii) dimuon + 2jets+ /Et, provided that the 
jets/leptons satisfy the requisite cuts. 



We present some numerical estimates where the five degenerate squark flavours are as- 
sumed to have masses ~ 400 GeV. The mass of the lighter stop is made to vary beween 
80 and 150 GeV, while the SU(2) gaugino mass M2 is held fixed at 150 GeV. Unification 
of the gaugino masses has been assumed. In addition, tan (3, the ratio of the two Higgs 
expectation values, has been fixed at 3. Although we have not assumed any definite high 
scale SUSY breaking mechanism, it can be verified that such combinations of parameters 
can indeed be realised in a supergravity (SUSY) framework. This is possible by using one's 



freedom [JT3 ] with the trilinear soft SUSY breaking term A (while still preserving charge and 



colour invariance |14] ) and using the Higgsino parameter \x as a phenomenological input (|i~5 



something that is justfied if the Higgs mass parameters retain the freedom of differing from 
the 'universal' sfermion mass parameter (m ) at the grand unification scale. 

The above specified set of parameters allows one to calculate the R-conserving stop 
decay width. For the R-parity violating two-body decays, we have used three sets of values 
for A233 = A333, consistent with the expectation from the SK results. Using these, it is 
straightforward to calculate the branching ratios of the various stop decay modes. It is 
found that for the region of parameter space under investigation here, the branching ratio 
for R-parity violating decays ranges from 75 to 99 per cent, with a near-equal share between 
the br and b\i channels. 

A parton level Monte Carlo calculation has been performed for both single muon + 
three jets+ /E^ and dimuon + two jets+ /E^ final states arising from stop pair-production 
at the Tevatron. In the former case, we have demanded that at least one jet be identified as 
a ^-induced one. Both the top- and stop-production cross-sections have been QCD corrected 



using a a so-called i^-factor of 1.4 in case of top [|16j and 1.3 []17j for the stop. For hadronic 
decays of the tau, we have considered only modes with one charged track, arising from 
r^- — > 7t =1= 7Vt-, p ± u T , afu T . In order to calibrate our parton- level results, we have computed 
the numbers of the same types of expected signal events from top quark pair production. 
The numbers of both single-muon and dimuon events thus calculated in the latter case agree, 



within small errors, with the actual observations [y|,[n|] , so that our results may be treated 



as reasonable estimates of how many stop-induced events can be contained in the top signal. 

Since our main purpose is to see how stop signal can percolate into top signals for which 
detectors are already designed, both have been subjected to the same set of cuts, as specified 
by the CDF top search strategy In both types of final states of our interest, jets are defined 
using the cone algorithm as discussed in Ref. [ 375 ,|TD | . Jets are counted in the analysis only 



if \r)\j e t < 2 where r\ is the pseudorapidity |18|JT9|] . For single muon final state we require 



jets with Et > 15 GeV, at least one of which is identified as a 6-jet (with identification 



efficiency »s 40 % QlS | ) while for dimuon final state two jets with Et > 10 GeV are required 
|T9| . Two jets are merged if AR < 0.4, where AR = y/Arj 2 + Acf) 2 , Ai] and Acfi being the 
separations in pseudorapidity and azimuthal angle. In both the cases isolated muon(s) with 
Pt > 20 GeV are necessary in the central region with \r]\ < 1 ]n],[nj. The criterion for muon 



isolation imposed here is that the total Et within a cone of Ai? = 0.4 around it should be 
less than 10 % of its own pt- In dimuon final states, backgrounds from real Z are rejected 
by requiring the dimuon invariant mass to be outside the interval 75 — 105 GeV ||19||. For 



single muon events, missing transverse energy, /Et > 20 GeV is demanded [18] while for 



dimuons /Et > 25 GeV is generally required |1| . It is observed in the last of Ref. JT9| that 
this combination of cuts effectively eliminates backgrounds from Drell-Yan production and 
other relevant processes for the dimuon final state. 

The net contribution of the stop to the top-like signals depends on two effects: the 
production cross section of a stop pair vs that of a top pair, and the relatively large branching 
fraction for the stops cascading into the final states of our interest. In addition, of course, the 
susceptibility to cuts plays a role. For ~ 100 GeV the production rates approximately 
match, and the stop production rate falls for higher stop masses. 

Figures 1 and 2 contain our main results, based on an integrated luminosity of 109 pb _1 . 
We have checked by varying the parameter fi between —700 and +700 GeV that the results 
are not altered by more than about 25 per cent. The dimuon events enable us to rule out a 
slightly larger range of the stop mass. This is because of more suppression of the top signal 
in this channel through leptonic branching fractions of W as opposed to what happens for 

6 



stop-pairs. It is clear from figure 2 that the dimuon data should definitely rule out the 
entire stop mass interval that is otherwise constrained assuming B(ti — ► c%°) — 1. In fact, 
the limit can perhaps be pushed somewhat higher up in the R-parity violating case. For 
A233 — A333 = 10~ 4 , this limit should be around 125 GeV with the currently available data, 
while for A -values of 5 x 10~ 4 it touches about 140 GeV. 

From the single muon signal too, stop masses upto about 120 — 125 GeV seem to be ruled 
out. We also want to emphasize that both this channel and the dimuon one can be used 
effectively in the context of Run II of the Tevatron to look for SUSY signals burried within 
the top data. A careful analysis of the dimuon versus dielectron data there, invigorated 
by the additional available luminosity, may lead to successful identification of an excess in 
the former, thereby indicating R-parity violation of a kind that simultaneously explains the 
atmospheric neutrino puzzle. 

For simplicity, we have confined ourselves in the above discussion to cases where the stop 
is kinematically forbidden to decay into a b quark and a chargino. Such a decay can be 
possible when M 2 and consequently the mass of lightest neutralino is smaller. However, in 
view of the fact that the chargino can subsequently decay into a lepton, the above mode can 
in effect strengthen the stop contribution. 

The values of R-parity violating couplings which are as small as the ones under study 



here are unlikely to contribute to precision electroweak data [£0] or to R-parity violating 
decays of the top quark |21j such as t — > bl or t — > bl even if they are kinematically 
allowed. However, with couplings several times larger in magnitude (with larger average 
squark mass, for example), some of these decays could lead to additional signals (final states 
of higher multiplicity) from the decays of the top, particularly at Run II of the Fermilab 
Tevatron. 

Another signal, similar in nature to the ones discussed above, is one where both the 
stops decay into the b\i channel. In this case one shall see dimuons with two jets. Proper 
6-identifaction in the jets, together with effective measures to eliminate the Drell-Yan back- 
grounds, can establish such final states as rather effective ones for discovering a light stop 



in an R-parity violating theory. 

As has been mentioned at the beginning, the above discussion has assumed only trilinear 
R-parity violating interactions. The presence of the bilinear terms of the form L1H2 in 
the superpotential (and the consequent vacuum expectation values for sneutrinos) can give 
rise to mixing between charged leptons and charginos [22|, which can again lead to top-like 
signals of the same kinds from stop decays. A detailed analysis of the signals in such a case 



will soon be presented by us |23 



In conclusion, the first signals of R-parity violating SUSY observable in current and 
upcoming experiments can very well mimic those of the top quark. This is particularly true 
if the theory has to account for the neutrino masses and mixing suggested by experimental 
results on atmospheric muon neutrinos. We can use the already available single leptons jets 
as well as dilepton data from the Fermilab Tevatron to constrain a large range of the lighter 
stop mass in such a scenario. In the upgraded version of the Tevatron, too, top quark signals 
will continue to be useful probes for such types of SUSY. 

We thank S. Raychaudhuri and M. Guchait for technical help. 
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FIG. 1. Contributions to single muon + 3 jets + JEt events from tt production for mt = 175 
GeV (solid line) and t\t\ production (followed by R-parity violating decays of ti's, dashed lines) 
as a function of at Tevatron with y/s = 1.8 TeV and an integrated luminosity of 109 pb _1 
and using the same set of cuts (see text). The MSSM parameters used are : M2 = 150 GeV, 
fi = — 400 GeV, tan/3 = 3, # t ~(stop mixing angle) ps —45°, mQ^u,D = 400 GeV for first two 
generations of squarks. We take A' = A233 = A 333 . CTEQ-4M parton distributions are used. 
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FIG. 2. Same as in Fig. 1 except for contributions to di-muon + di-jet + JE? events. 
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